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Abstract. We calculate the probability distribution of the overlap between a 
spin glass and a copy of itself in which the bonds are randomly perturbed in 
varying degrees. The overlap distribution is shown to go to a i distribution in 
the thermodynamic limit for arbitrarily small perturbations (bond chaos) and we 
obtain the scaling behaviour of the distribution with system size TV in the high 
and low temperature phases and exactly at the critical temperature. The results 
are relevant for the free energy fluctuations in the Sherrington-Kirkpatrick model 

ffl- 



PACS numbers: 75.50.Lk, 75.10.Nr 
1. Introduction 

Chaos is one of the fascinating properties of spin glasses. Whenever a spin glass is 
subjected to a change (e.g. a change in temperature or in the coupling constants), the 
equilibrium state of the system changes completely (in the thermodynamic limit), no 
matter how small the change. This behaviour has been first observed in hierarchical 
models [2] and was then suggested for finite dimensional spin glasses |3] within the 
droplet theory. Chaos is however not restricted to replica symmetric droplet-like 
scenarios: in 4 it was shown that temperature chaos exists in the Sherrington- 
Kirkpatrick model [5] which breaks the replica symmetry. In this paper we show that 
bond chaos also exists in the Sherrington-Kirkpatrick model. This is not surprising 
since bond chaos is usually a stronger effect than temperature chaos [6]. However, it 
has recently been found that an intricate connection exists between bond chaos and 
the sample-to-sample fluctuations of the free energy [lj [7] , which is a long standing 
problem in spin glass and extreme value theory. It is therefore necessary to calculate 
bond chaos in order to make progress on fluctuations, and in this paper we will present 
the details of the calculation which was only sketched in pQ . 

Temperature chaos has been simulated extensively in the literature, see e.g. [5] and 
references therein. Bond chaos has not been simulated quite so much but nevertheless 
for many different models, including the Edwards- Anderson model in dimensions 1 
to 4 |1 EH ED U, on the hierarchical Berker lattice [U [T3] and on Bethe lattices 
[TO] but, to the best of our knowledge, never for the Sherrington-Kirkpatrick model. 
Here we will show that the probability distribution of the overlap, which is our main 
quantity of interest, has a relatively complicated finite size scaling behaviour (see 
Eq. d67]) below). 

This paper is organized as follows. In Sec. [2] we specify the model and method we 
will be using and explain in Sec. how to obtain the probability distribution of the 
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overlap for this model. We then proceed with the replica calculation in Sec. 0] where 
we explicitly calculate this distribution above, at and below the critical temperature. 
We end with a conclusion in Sec. 



2. Model 



We will compare two real replicas of a Sherrington-Kirkpatrick spin glass with 
Hamiltonian 

H e = -—=^2 K ii( e ) S i S j' C 1 ) 

i<j 

where the Sj (i = 1, . . . , N) are N Ising spins and 

are Gaussian random variables of unit variance, composed of independent Gaussian 
random variables Jij and J^ , also of unit variance. The parameter e is a measure of 
"distance" between the sets of bonds {Kij(0)} and {.Ky(e)}. If e = 0, the bonds are 
equal, if e — oo, the bonds are completely uncorrelated. 

Our main question concerns the disorder averaged probability distribution of the 
spin-spin overlap P e (q) between the two replicas, the first with Ti.o and the other with 
H t . Here q is defined by 

1 1,0 2,e /on 

i 

where s^' x is the ith spin in replica r, which has Hamiltonian TL X . 

We know that for temperature chaos (i.e. when comparing two replicas with 
identical bonds but at different temperatures), even an infinitesimal temperature 
difference AT leads to Pat{q) — $(q) in the thermodynamic limit [3]. This means 
that equilibrium states in the two replicas are totally uncorrelated. Since bond chaos 
is usually a stronger effect than temperature chaos, we expect the same here, and we 
will show this explicitly below. More interesting, however, and also more important, 
is the question of how P e (q) scales towards the 5-function with system size N. This 
information is for instance needed for the calculation of the sample-to-sample free 
energy fluctuations [H [7] . 



3. Probability distribution of the overlap 

In order to calculate P e (q) we first try to formally calculate P e> j(q), the nonaveraged 
probability distribution to find the overlap q. Here and in the following, the subscript 
J indicates a nonaveraged quantity. Given a realization of the disorder for the two 
replicas, we can write down a partition function Z e ^j(q) for them, constrained to have 
the overlap q, 

Ze,M = Tr5 (q- ^E S ''°V J ex P (/?£i^(0)^ 4'° + /3£^)s-' e 4 £ I -(4) 
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This method of constraining two systems to have a given overlap was first suggested 
in [2]. From this one gets the free energy j3F t j{q) = —logZ e j(q), and defining 
Y e> j = f dq Z e ^j(q), the probability distribution of q is given by a Boltzmann factor, 

p-PF Ci .,(q) z ,(a) 
P eJ (q) = = (5) 

Of course, we cannot calculate F €t j(q) for a given disorder but we can calculate its 
disorder average F e (q) — E F e> j (the disorder average is denoted by the symbol E ■ ■ ■) 
in the thermodynamic limit by replica methods. This will be done below. But first 
we want to continue formally in order to assess the approximations we are going to 
make. 

We split F e> j(q) into three parts, 

Fe,j(q) = Nf e (q) + AF,(q) + AF eJ (q). (6) 

The first part, Nf e (q), is the extensive part of the average free energy. The second 
part, AF e (q), is the finite size correction to it. Finally AF ej j(g) is the fluctuation of 
the free energy about its disorder average F e (q), so E AF ey j(q) = 0. 

In order to calculate P e (q) = E P e j(q), we first look at Y e j. It can be written as 
follows, 

Ye j = f 1 dqe -PNf e (q)-f3(AF e ( q )+AF e: j(q)) 



(-/?)" ^dqe-^ N f^\AF e (q) + AF e „,(g))" \ 
n\ jj dqe-^l-M) j 

= Y ? ( 1 + E t(AF e + AF e>J )"] ) , (7) 



nl 

71=1 



where Y" e ° := J^dqe & N f^i) and [• • -]o is the average taken with the probability 
distribution 

p-/3JV/ e (g) 

P °M = — po— ( 8 ) 



We then have 



and 
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P t (q) =P°{q)E f&e—E: 1 ^[(A Fl+ AF ( . J )-] -^ e (,) + AF ( , J (,)) (1Q) 

Jo 

= P°(g)(l-/3(AF £ (g)-[AF £ ] ) 

[(AF e (g) - [AF^)\ - ^E [(AF e>J (q) - [A^, 7 ] ) 2 ] 

+ ^ (AF e (q) - [AF t } Q f + ^E (AF tJ (q) - [AF e J Q ) 2 + •■■). (11) 

The last line follows from expanding the exponentials and sorting by powers of AF e 
and AF £jJ . 

Eq. (JTTJ) shows precisely what can be expected of the calculation which is to 
follow. We will calculate the first term of this expression, P®(q). This is the 
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same term one would have written down from the start by following large deviation 
statistics principles. Anything beyond this would require us to calculate the finite 
size corrections AF e (q) and AF e ^(q) which is impossible due to the massless modes 
present in the spin glass phase. But our precise analysis allows us to gauge the 
applicability of the large deviation statistics approximation: although AF £ (q) and 
AF e j(q) are not necessarily small (they grow with some subdominant power of N), 
they only appear in the combinations AF e (q) — [AF e ]o and AF €} j(q) — [AF fj j]o. If 
these differences are small, our results not only hold for "large" deviations but also 
for "small" ones. While we have no proof for this, application of our results to the 
free energy fluctuations above and at the critical temperature [1] [7] shows that it is 
at least true at those temperatures. As to the low temperature phase, we will see 
later that for e = and q less than the Edwards Anderson order parameter ^ea, the 
corrections are of order 1 but small enough not to introduce any qualitative changes, 
and we expect the same to hold for e > 0. 

4. Replica calculation 

Temperature chaos in mean- field spin glasses has been treated in the literature [15j [4] 
and we refer the reader to these papers for details. Repeating Rizzo's calculation [15] 
(see also [TB] for the calculation at e = 0) but for bond chaos rather than temperature 
chaos, one arrives at the following expression for the disorder averaged, replicated, 
constrained partition function E Z™ j(q) 



Here, r = ~(1 - 1//3 2 ), t' = ±(1 - \/l + e 2 / ' (3 2 ) and the 2n x 2n matrix T is given by 



and the n x n matrix Q is zero on the diagonal, while R, also of size n x n, is not. 
The latter will therefore be split into a part on the diagonal pdX and the rest P. At 
this stage it is an ansatz to make all entries on the diagonal of R equal to pa but this 
need not be done and one could proceed without this assumption and justify it later. 

Note that, in contrast to [15], we are using the truncated model here by keeping 
only those fourth order terms which are written down in Eq. (|12p . This will make the 
calculations below much easier. 

The saddle point equations following from this expression can easily be derived 
and, making a Parisi symmetry breaking ansatz for P and Q, one gets the following 
system of equations: 




(12) 



) 




(13) 



= rq(x) + V(.*) f"dz (q(x) - q{z)f - ~ f* dz (p(x) - p(z)) 2 




= r'p(x) +-p (x)-w I dz (p(x) - p(z))(q(x) - q(z)) 
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-w(p ~ p d )p{x) ~ wqq(x) (14) 

-~w(p - Pd)q{x) - wqp(x) (15) 

2y f 1 

q = Pd + YPd- 2w J dzp(z)q{z). (16) 

The variable q (without argument) in the last equation is the overlap the two real 
replicas are forced to have. As usual, p and q denote the integrals over p{x) and q(x) 
from to 1. The free energy per spin may be expressed in terms of the solutions p(x), 
q(x) and pd of these equations and is given by 

M ( q )= qPd -§-f- + rj\z q 2 (z) + lj\z q\z) 



dzzq 3 (z)—w / dzq 2 (z) / dz' q(z') + 2wpd / dzp(z)q(z) 



w 

+t' [ dzp 2 (z) + - I dzp A (z)-w I dzp 2 (z) I dz'q{z') 



It can be checked that this free energy gives back the saddle point equations Eqs. (114 
(|16|) when the derivatives with respect to p(x), q(x) and p d are taken. 



4--1- Above the critical temperature 

Above the critical temperature, r < 0, the saddle point equations can be solved 
perturbatively in the limit of small q. For q = 0, the exact solution is q(x) — p(x) = 
and P d = 0. For < q <C |t| it is easy to see that q{x) — 0(q 2 ), p{x) = 0(q 2 ) and 
Pd = q + 0(q 3 ). Plugging this into the free energy, Eq. (fT7|) . yields 



2 V VI 

The probability distribution P®(q) is thus 

F e °(g) oc e ~N(4h(e) + o(^)) (19) 
with h(e) = 1 - -JL-. 



4-2. At the critical temperature 

At the critical temperature, r — 0, we can solve the saddle point equations in the 
limit q -C \t'\. Just like above the critical temperature, the solution is q(x) = 0(q 2 ), 
p(x) — 0{q 2 ) and pd = q + 0(q 3 ), and the free energy is also given by Eq. (fT5| in this 
limit. Expanding in powers of e, we find 

Pfe(q) = jh(e) = ^f+ O(e^) (q « e 2 « 1). (20) 

We will also need the solution of the saddle point equations in a different limit, namely 
e 2 < 5 « 1. To obtain this, it is convenient to rewrite the saddle point equations and 
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the free energy in terms of the functions a(x) := q(x) + p(x) and b(x) :— q(x) — p(x). 
A straightforward calculation leads to 

= (r + wp d )a{x) - T -^(a{x) - b{x)) + j^(a 3 (x) + 3a(x)b 2 (x)) 



dz (a(x) — a(z)) 2 — waa(x) (21) 



= (r - wp d )b(x) + l—J-{a{x) - b(x)) + ^(b 3 (x) + 3b(x)a 2 (x)) 

_^ f dz [b(x) ~ b{z)f - wbb(x) (22) 
2 Jo 

9=W + Y^-f J dz(a 2 (z)-b 2 (zj) (23) 
and 

PMq) = QPd-§-^- 2(1 9 -2t<) + lj! dz (fl2(z) + h2[z)) ~ T ~T~ t dz {a{z) ~ b{z))2 

-- [ dzz(a 3 (z) + b 3 (z))-- [ dz [ dz' (a 2 (z)a{z') +b 2 (z)b(z')) 
6 Jo 2 Jo Jz 

/ dz(a 4 {z) + 6a 2 (z)b 2 (z) + b 4 {z)) + ^. f dz (a 2 (z) - b 2 (z)) (24) 
48 Jo 2 Jo 

These saddle point equations can easily be solved for r = t' = 0, and the solution is 



( 2wx x < x' 

a{x) = \ 2 L' 2 . ? (25) 

- X > Xo 



b{x) = (26) 
2y 3 2w\ , 2 2 , 3 , 



^ 3 eon , ,2 * ,3s 

<7 = Pd + — P d - — a" (a; 2 - -x 2 ). (27) 



The breakpoint of a(x) is x' 2 = 1 — yl — The value of has to be obtained by 
solving Eq. (f27|) for p^. 

When we want to find the free energy perturbatively in the limit < |r'| <C g <C 1, 
we notice that Aa(x) and Ab(x) (the correction terms for nonzero r') are both of order 
t'. Since these are the corrections to the saddle point solution found at t' — 0, they 
only contribute to the free energy to order r' 2 and can therefore be neglected as long 
as we are only interested in the free energy to order r'. Plugging the functions a(x) 
and b(x) just found into Eq. (124]) yields 

0fe(q) = ^q 3 --r'q 2 +0(r' 2 )+0( q 4 ) = ^+^ q 2 +0(e\ q 4 ) [e 2 « q « 1).(28) 
6 4 6 16 

We notice that for e <C iV -1 / 6 the exponent Nf3f e (q) of the probability 

distribution P e (q) ~ e - N 0M<i) [ s dominated by Nwq 3 /6 from Eq. ([28]) . while for 

jy-i/6 < e < i it is dominated by N\h(e) from Eq. ([20]). Therefore the probability 

distribution is 

P M S e- Nw « 3 / e e « AT-V6 
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4-3. Below the critical temperature 

4-3.1. Some exact solutions of the saddle point equations Eqs. (|14p - (|16[) can not be 
solved for general q and e (or rather r' which contains the only reference to e). They 
can however be solved exactly for e = and for q = and perturbatively in various 
limits. 

The solution for q = is simply 

Pi(as) = (30) 

qi(x) = \ X T>r (31) 

Pd =0 (32) 



where x 2 = 1 - J 1 - |£ . That this is a solution can easily be checked since Eq. (fH)]) 

is clearly satisfied, Eq. (fT5|) is also satisfied by pd = and p(x) = 0, irrespective of 
q(x), and Eq. (TT?]) in this case reduces to the normal Parisi equation the solution of 
which is well known and given in Eq. (|32[) . 

The solution for r = r' and q < ^ea, where <7ea is the Edwards- Anderson order 
parameter, has been found by Rizzo [15] and is given by 

J" wx x < Xl 

P2(x) = { ^ x > Xl ( 33 ) 



II 

wx 
V 

If. I ] 

y 
2y 

It' J- 2 

2y 



x < Xl 

xi < x < 2xi , 

o / 34 

LXi < X < X2 
X2 < X < 1 



Pd = T ^- r (35) 

with xi — ^f-. This solution however only exists if 2xi < X2, which is the case when 
q < 9ea- In the limit q — > it goes to the solution found above for q = 0. As shown 
by Rizzo, the free energy /3/o of this solution (for any q where this solution exists) 
is precisely equal to that of two unconstrained systems; we may therefore use /3/o as 
reference energy. 

For r = t' and q > (?ea, another solution takes over [16j . It is given by 



LUX 



Ps(x) = q 3 (x) = { ^ ^ > ^ ( 36 ) 



y 

2y 3 2w 3 2 2 c 



Pd + ^-P d -—r(xi--xi). (37) 



The breakpoint is x 3 = 1 — y 1 — yT+ ™? Pd , and can be obtained by solving Eq. ([37]) 

for Since the details can be found in [16] , we merely quote here that the excess 
free energy of this solution is 

0f-0fo =c (q-qEA) 3 (38) 
with some positive constant cq. 
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4-3.2. Perturbative solution in the limit e — > oo The first limiting case to consider 
is e — > oo or, equivalently, t' — > — oo. In this case, dividing Eq. (| 1 5|) by r' yields to 
leading order 

wp d q(x) 



p(x) 



+ 0(1/t u ). 



Inserting this into Eq. (I14|) gives, again to leading order, 

= ( T + ^)q(x) - V -q\x) -^fdz (q(x) - q{z)) 2 



(39) 



wqq(x),(AO) 



i.e. precisely the Parisi equation but for a slightly lower temperature r + w _^f ■ 

Let's look at the free energy of this solution in the limit r' = — oo. From Eq. (UHJ) 
we get qpd = p 2 d + 2yp\jZ and the free energy difference /3/oo(<z) — flfo 1S 

(41) 



yp d 



This means that the probability density Poo(q) is proportional to exp(— N(q 2 /2 + 
C(g 4 ))), from which [q 2 ]o can be evaluated and is given by [q 2 ]o = j^. This is a useful 
result to check whether the method works since it is easy to convice oneself that this 
is indeed true for two replicas with completely uncorrelated couplings: 



E(q 2 ), 

in agreement with the replica result 




. i<3 



4-3.3. Perturbative solution for e 2 <C q <C 1 We now turn to the opposite limit, 
6 — ► or t ! — ► t. We first consider the case e 2 q 1. In order to solve the saddle 
point equations in this limit, it is again convenient to consider the version in terms of 
a(x) and b(x), Eqs. (f2Tj) - (j23j) . The solution for r' = r, i.e. At = e = 0, and arbitrary 
q (small enough such that the solution exists), taken from Eqs. (|35p - (|35|) , is 

q 



a(x) 



b(x) 



- 2t 

2wx 

y 

2wx\ 

w{^x+2x i) 

W(X2 + 2X! ) 



w{x — 2x±) 
2y 

w(x 2 — 2xj) 

2y 



X < X\ 

x\ < x < 2xi 

2xi < X < X2 
X2 < X < 1 

x < 2x\ 

2x\ < X < X2 
X2 < X < 1 



(43) 



(44) 



(45) 



Surprisingly, a perturbative expansion of the saddle point equations for small At shows 
that the leading order correction to this solution is of order \J At. It is straightforward 
to show that it is given by 

x < x\ 

AM.ri = { Xl <x<2 Xl . (46) 

2xi < x < 1 
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Inserting this into the free energy difference leads to 

Pm-Ph = if A, = 12 J^ 2t)3 =: Cl eV « a « 1), (47) 

defining the constant c\. 

4-3.4- Perturbative solution for q <C e 2 <C 1 In this limit we can use the solution 
for q = from Eq. (|32p as reference solution and construct a perturbative solution 
for small q (or p e i) from it, i.e. we set p{x) = pi(x) + Ap(x) — Ap(x) and 
g(x) = qi(x) + Aq(x) and write down the lowest order in pd from Eq. (fTS"]) : 



= r'Ap(a;)-«; / (Ap(x)~Ap(z))(q 1 (x)-q 1 (z))-w(Ap-p d )q 1 {x)~wqTAp{x).{A8) 
Jo 

As we can see, Aq(x) does not appear in this equation at all, and since p{x) only 
appears quadratically or in combination with p d in Eq. (TH)) , we conclude that Aq{x) 
must be in fact of higher order (0(p d )) such that we may set it equal to zero. 

Eq. (|4"5)) can be solved in the following way. First note that qi — t/w such that 
with At := r — t' 



= -AtA p (x)~w dz(Ap(x)-Ap(z))(qi(x)-qi(z))-w(Ap-p d )qi(x). (49) 
Jo 

Differentiating once with respect to x yields 

C-X C-X 

= -ArAp'(x) - wq[(x) / dz (Ap(x) - Ap(z)) - wAp'(x) / dz (qi(x) - qi(z)) 

Jo Jo 

-wt&p'-pMx). (50) 

We conclude that Ap'(x) = whenever q[(x) — 0. When q[(x) ^ 0, differentiating 
Eq. (fST))) once more gives (noting that g"(x) = 0) 

Q = -ArAp"(x) -2wxq' l {x)Ap'{x) -wAp"{x) [ dz {q x {x) - q x {z)). (51) 



9 2 \ 2 

At + J Ap"(x) —Ap'(x). (52) 

This differential equation can easily be solved and one gets 

Ap'(x)=C'(at + ^^ 2 (53) 
with an as yet undetermined constant C. Integrating once with respect to x yields 

^>-^'te)" c 'te)' (54 » 

where 

/(z) = t + arctanz. (55) 

1 + z 2 y ' 

This is true for x < X2 while for x > X2, Ap(x) stays constant. We can now calculate 

Ap, 

A~p = C r dxf(-^=)+C(l-x 2 )f(-^=) (56) 



V2 v ^at; 



= C arctan(^ 2 /2 V^A7) + (7(1 - z 2 ) WX% ^^f\ ■ (57) 

1 + w z a;2/4yAr 
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The constant C can be determined from Eq. (|50|) by taking the limit x — > 0. In 
this limit the equation reads 

ArAp'(O) = w(p d - 3£)gJ(0) (58) 

such that 

Cwv/Ar w 2 w 2 -r— 

which yields upon expanding in powers of V At, 

C=VfW|- *») f + ^(Ar 5 / 2 ). (60) 

7T 7T^ \3 / \ WX2 J 

With this value of C, p(x) = Ap(x) together with q(x) = q%(x) is the correct solution 
of Eqs. p4p and (fT5|) for small Ar up to order pd- 

What of the free energy of this solution? In order to calculate this, we eliminate 
q from Eq. (fl~7|) using Eq. (fT6)l and subtract /3/o to get 

Pfe(q) - Afo = f + ^ - 2{1 ^ 2t , } (Pd + ^f- 2w£ dzp(z) qi (z) 

-At f dzp 2 (z) + ^ f dzp 4 (z)-w [ dzp 2 (z) f dz' {q x {z) - q x [z')) 
Jo 6 J Jq J 

-2w [ dzp(z) qi (z) [ dz'p(z'). (61) 



We can calculate the integrals appearing in this equation and obtain 

A f \ I \ Tpd Arpd i ^Vv AT3/2 Pd ,, nfA 5/2n , R0 ^ 
dzp(z)qi(zj = 1 t. h C(At j/ ) (62) 

W W TTW 

rl dzp 2 (z) =P l-^lpl^Ar + 0(Ar^) (63) 
Jo ™ 
and 



dzp(z)U(z) J dz' (q 1 (z)-q 1 (z')) + 2q 1 (z) dz' p{z')^j 



= L A _ 2 ^ A + ^S^!l A + o ( ^). (64) 
w w irw z 

Assembling these results, the terms of order Ar cancel and the free energy difference 

is 



Pfe(q)~Pfo 

TTW 

_ A^y 1 1 
C2 -^(l-2r) 2 W { ) 

4-4- Perturbative solution for q -C min(l, e ) 

In the preceding subsection we have calculated the free energy for 5 < e < 1. The 
condition e <C 1 is an unnecessary restriction, however, and one can in principle 
carry out the same calculation as above without expanding for small e, as long as 
q <C min(l, e 2 ). The result is 

PMq)-Ph = f(t)q 2 + 0(q 3 ), (66) 
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where the function /(e) is monotonic and has the properties /(e) = c^e 3 + 0(e 5 ) for 
e — > and /(e) — ► | for e — > oo. The latter follows from the solution for t' — > — oo 
we found previously. As we won't need any more detailed information about /(e), we 
will not carry out this tedious calculation in detail here. 

4-5. Probability distribution 

As before, we want to calculate the probability distribution P®(q) ~ e ~ N P(f--X<i)-fa) _ 
We observe that the admissible range of e (the interval from to oo) divides into 
4 parts. For e <C TV -1 / 2 , both Eq. (|4"T|) and (p5|) produce a negligible exponent 
N (3(f e (q) — /o) f° r all q G [0, <?ea]- The probability distribution P e (q) is thus a constant 
in that interval, with an exponentially suppressed tail for q > <;ea from Eq. (|38| . In 
the range iV" 1 / 2 <C e <C TV -1 / 5 the exponent N[3(f <L (q) — / ) is negligible for q <C e 2 
when Eq. |65|) prevails. It only becomes noticable when q > N~ 2 ^ 5 e 2 such that we 
can approximate it by Eq. (|47p . For N^ 1 / 5 <C e < eo with some arbitrary small and 
fixed eo independent of N, the exponent is dominated by Eq. (|65|) . Finally, for eo < e 
we can use Eq. (|66|) . 

Combined, we can write 

% - «ea) e « JV-V2 

e-RTc^V jy-i/a < £ #-1/6 

e-^«V 7V- 1 /5« e < eo ' (b7) 

e-^M? 2 e < e 

where the function 9{x) = 1 for x < and #(x) = e~ Nc ° x3 for x > 0. In principle, 
the latter two regimes could be combined to P £ (q) oc e^ N ^ e ^ q for N^ 1 ^ 5 <C e but 
splitting them makes the dependence on e more explicit. 

A note is in order about the probability distribution at e = 0. Clearly, the result 
in Eq. (|67[) does not coincide with the known distribution from the Parisi solution 
[T71 ri"8l IT9] , in particular the 5 peak at (?ea is missing and in the range < q < qea the 
probability distribution is flat. The latter is not only an artifact of using the truncated 
model (for which the distribution would indeed be flat) but was already derived in 
[To} ITB] for the full model. Both the flatness and the missing S peak originate from 
neglecting the finite size corrections in Eq. (jlip. We conclude that these corrections 
are of order 1 for q < qea (but do not change the probability distribution qualitatively, 
i.e. do not introduce gaps or zeros in the distribution) and conspire to form the S peak 
for large N at q = qea- When e > N^ 1 / 2 we do not expect any 5 peaks since the 
equilibrium states in the two replicas start to differ substantially, and in the light 
of what we saw at e = we expect the finite size corrections to be similarly good- 
natured. Thus the results we have presented here for the probability distribution 
should be correct up to prefactors which might vary for the full model. 

5. Conclusion 

We have shown that bond chaos exists in the Sherrington-Kirkpatrick model by 
calculating the probability distribution of the overlap q between two copies of a system, 
one of which has randomly perturbed bonds with respect to the other. The finite 
size scaling of this distribution has been calculated above, at and below the critical 
temperature. In the low temperature phase four different regimes have been identified. 
For bond distances e <C iV -1 / 2 the distribution has a variance of order 1, i.e. the 
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equilibrium states in the two copies are still very similar. For TV x l 2 -C e -C N 1 / 5 , 
the distribution is proportional to e~ Ncit q , i.e. already very narrow, the width scaling 
as ^ 2 ^ 3 with the scaling variable £1 = \f~Ne. For TV -1 / 5 <Ce<l, the scaling variable 
is £ 2 = TV 1//3 e and the distribution becomes Gaussian with width proportional to £^~ 3 ^ 2 . 
For all other values of e, finally, the distribution remains Gaussian, and its width goes 
as JV-V2. 
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